Short-ranged RVB physics, quantum dimer models and Ising gauge theories 
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Quantum dimer models are believed to capture the essential physics of antiferromagnetic phases 
dominated by short-ranged valence bond configurations. We show that these models arise as par- 
ticular limits of Ising {Z2) gauge theories, but that in these limits the system develops a larger local 
f/(l) invariance that has difi'erent consequences on different lattices. Conversely, we note that the 
standard Z2 gauge theory is a generalised quantum dimer model, in which the particular relaxation 
of the hardcore constraint for the dimers breaks the (7(1) down to Z2. These mappings indicate 
that at least one realization of the Senthil-Fisher proposal for fractionalization is exactly the short 
ranged resonating valence bond (RVB) scenario of Anderson and of Kivelson, Rokhsar and Sethna. 
They also suggest that other realizations will require the identification of a local low energy, Ising 
link variable and a natural constraint. We also discuss the notion of topological order in Z2 gauge 
theories and its connection to earlier ideas in RVB theory. We note that this notion is not central to 
the experiment proposed by Senthil and Fisher to detect vortices in the conjectured Z2 gauge field. 



PACS numbers; 75.10.Jm, 74.2Q.Mn 
I. INTRODUCTION 

The question posed by high-temperature superconduc- 
tivity is how |a Mott insulator becomes superconducting 
upon doping. EJ As the insulator is itself, at low energies, 
also an antiferromagnet hostile to the motion of holes, 
much work has been based on the notion that the doped 
state is best approached from a "nearby" insulating state 
that lacks long range order, i.e. a spin liquid.El The dopejd 
spin liquid is then argued to become superconductingj3 

The simplest such scenario casts the resonatingpvalence 
bond (RVB) state proposed in 1973 by AndersonP in the 
role of the spin liquid. Pairs of electrons form singlet 
(valence) bonds, a superposition of which yields a liq- 
uidlike, non-Neel ground state. Holes doped into this 
state undergo spin-charge separation. The charge degrees 
of freedom, able to move freely through the spin-liquid, 
become superconducting upon Bose condensation. The 
spin excitations are understood as composites of spin- 
i spinous and the decay of the electron into holon and 
spinon provides a natural explanation of the broad quasi- 
particle spectra seen over much of the normal state of the 
cuprates. 

RVB scenarios themselves cover a broad range of possi- 
bilities. The short-ranged (SR) flavour of RVB stays close 
to Anderson's original vision by including valence bonds 
only between electrons located in a small neighbourhood 
of one another leads to gapped spinonsQ. Its low en- 
ergy dynamics is believed to be captured by the quan- 
tum dimer model (QDM) introduced in Rcf. ||, where a 
VB is represented by a dimer linking the two electrons 
which form it. Historically, the short-ranged RVB was 
abandoned when the QDM failed to lead to a spin liquid 
on the square lattice - it typically leads to a columnar 
state - and was considered suspect for building in a spin 
gap (i.e. a gap to triplet excitations) that was not in evi- 
dence at optimal doping; subsequent to the identification 



of the pseudogap regime and the discovery of stripes these 
"defects" seem less compelling although the problem of 
describing the collapse of the spin gap with doping is still 
unsolved in this approach as is the still basic probleM of 
solving for the physics of a finite density of dopants. Q 

In contrast the long-ranged RVB versions have received 
gauge theoretic treatments based on U{1) and SU{2) 
reformulations of the Heisenberg model tha.t can give 
rise to a gapless mean-field spinon spectrum.cl While the 
broad similarity between the mean-field phase diagrams 
constructed early on, and the phase diagram of generic 
cuprate superconductors is striking, assessing the impact 
of fiuctuations has been difficult. In particular, the gen- 
eral belief that such gauge theories cannot give rise to 
deconfined phases in 2-1-1 dimensions is at odds with the 
program of finding a proximate fractionalized spin liquid. 

Recently, Senthil and, Fisher (SF),El building on earlier 
work by Balents et a/. ,113 have proposed to get around this, 
by reformulating the problem as an Ising gauge theory.Ell 
As Ising gauge-Lheories do have deconfined phases in 2+1 
dimensions this seems quite promising. What is not 
clear from their work, is exactly what microscopic degrees 
of freedom are described by the Ising gauge fields .O SF 
have offered two separate justifications for the presence of 
Ising gauge fields. First, that a four fermion Heisenberg 
interaction can be decoupled by means of an Ising gauge 
field and second, that in models with separate electronic 
and superconducting degrees of freedom, the latter can 
screen the charge of the former upto a sign ambiguity in 
defining the needed square root of the cooper pair oper- 
ator. The former seems to us to be an interesting and 
exact microscopic statement, but inconclusive regarding 
the nature of the low energy theory; this point has also 
been made recently by Hastingst£l and we will give a triv- 
ial example to illustrate this point later in the paper. 
The second justification builds in the physics invoked 
in earlier work, namely the capacity of a superconduct- 
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ing condensate to screen charge and turn quasiparticles 
into spinonsJlj but it does appear puzzling that it holds 
into insulating phases as hypothesized by SF. In addi- 
tion SF have argued that a deconfined phase involving 
Ising gauge fields must be characterized by the notion of 
topological order invoked in studies of the quantum Hall 
effectim and that this order can be directly detected in 
an experiment. 

In this paper we attempt to shed some light on these 
issues by showing that at least one realization of Ising 
gauge physics is exactly the physics of the short ranged 
RVB. We will do this by formulating the QDM descrip- 
tion of the latter as an odd Ising gauge theory, a term 
we will explain in Sect. [II below. In this we will of- 



fer a variation on previous work by Fradkin and Kivel- 
son whojpapped the problem onto a [/(I) gauge theory 
instead.EEl This builds on a completely differently moti- 
vated stream of work of two of UR-Lwith P. Chandra) on 
frustrated quantum Ising modclst£l which has included 
a recent demonstration that the quantum dimcc_model 
on a triangular lattice supports an RVB phase.Ej This 
connection will allow us to interpret various statements 
about Ising gauge theories in the language of valence 
bonds - it will turn out that the Ising variable is the 
number of valence bonds - and, we hope, make them 
easier to grasp and evaluate. We should note that al- 
ternative identifications of spin liquid physics with Ising 
gauge theories in different limits have beea^-pade previ- 
ously implicitly by Read and Chakraborty,EM.a|iiil explic^ 
itly by Read, Sachdev, lalbert and Vojta,LJO Wen,E£l 
and Mudry and Fradkin .E3 

A second benefit of this exploration is that it focusses 
attention on what it takes to get an Ising gauge descrip- 
tion of the low energy dynamics, namely a binary link 
variable and a local constraint. If SF are correct and 
the Ising description has general applicablity, it should 
be possible to make comparable identifications in other 
contexts. 

In the balance of the paper, we will review the QDM 
description of valence bond phases, describe the refor- 
mulation of the QDM as an Ising gauge theory and of 
general Ising gauge theories as generalised dimer models 
(GDMs), collate the known results on these models, dis- 
cuss the notion of topological order in their context and 
conclude with a brief summary. As much of the interest 
of the paper lies in the connection between QDMs and 
Ising gauge theories, we have felt it useful to review a 
number of known results on both. It is useful perhaps 
to list the results that are new to this paper. These 
are the formulation of the QDM as an Ising gauge the- 
ory (Section III), the introduction of the odd Ising gauge 
theory and its QDM limit (Section IV) and its identifi- 
cation with the action including a Polyakov loop (PL) 
term written down by SF (Section IVB and Appendix 
B), the (elementary) treatment of Ising gauge theories in 
d — \ including the presence of fractionalization in the 
odd case, the rationalisation of the QDM results on the 
square and triangular lattices via the absence or presence 



of charge-2 Higgs fields (Section VB) and the discussion 
of topological order in Ising gauge theories (Section VI). 



II. VALENCE BOND PHASES AND QUANTUM 
DIMER MODELS 

Consider an insulating magnet with enhanced quan- 
tum fluctuations (as is the case with S" = 1/2 and other 
"small" values of the spin) and competing interactions 
that frustrate any Neel ordering that would be deduced 
from a semiclassical analysis. In the extreme case where 
any residual order is extremely short ranged, there is a 
reasonable expectation that the system will construct its 
ground state from configurations in which all spins are 
paired in nearest neighbour valence bonds. By continu- 
ity, we expect that there will be nearby Hamiltonians for 
which valence bonds of finite length will suffice and these 
are expected to share their basic physics with the purely 
nearest neighbour case.tZl As there is still a large number 
of short ranged valence bond states possible, even with 
the restriction to this sector there is a non-trivial problem 
remaining - that of diagonalizing the Hamiltonian within 
this highly degenerate manifold - which is the problem of 
"resonance" .□ Depending on the details of the Hamilto- 
nian, several phases might be realized. This set of phases 
are what we call (short ranged) valence bond dominated 
phases and, by hypothesis, they are all characterized by 
a spin gap. 

There are two primary obstacles to investigating the 
physics of valence bond dominated phases which we will 
restrict, in the remaining, to those with purely nearest 
neighbor bonds. The first is basic, namely the large de- 
generacy cited before, e.g. on the square and triangular 
lattices in d = 2 there are e°'^ (a > 0) states on an N site 
lattice. The second is technical in that different valence 
bond configurations are not orthogonal, although their 
overlap is effectively exponentially small in the number 
of dimers in which they differ. In some cases .-there is 
also a proof that they are linearly independent O 

To deal with the second problem it is convenient 
to formulate an expansion that can include the non- 
orthogonality perturbatively. As the parent configura- 
tions arc in one-to-one correspondence with hard core 
dimer coverings of the various lattices, such a tack leads 
to a quantum dimer model. 

The sites the electrons reside on define the direct lat- 
tice. The Hilbert space of the QDM thus consists of all 
hardcore dimer coverings of the direct lattice. The QDM 
Hamiltonian for the insulating case (half-filling) consists 
of two parts, a kinetic (T) and a potential {V) one. 
The former is off-diagonal and generates the resonance 
plaquette moves between different dimer configurations, 
whereas the latter, diagonal one, counts the number of 
plaquettes able to participate in such resonance moves. 

For the square lattice we find the Hamiltoniantl 



QDM 



= -tT + vV 
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En {-t 'I + h.c) + V + |I I) (I I|) } , 

(1) 

where we have kept only the simplest kinetic and poten- 
tial energy terms with coefficients t and v, and the sum 
Ed runs over all plaquettes. In what follows, we refer 
to the QDM with t; 7^ as the extended QDM, whereas 
QDM on its own refers to the case v — 0. On other lat- 
tices T will take the form of a sum of resonance moves on 
the shortest even loop (which is a plaquette in this case) 
and V will be count the number of such possible moves 
in a given configurations, e.g. on the honeycomb lattice 
both will involve three dimers. 

There are two main paths for obtaining a quantum 
dimer model from a magnetic system. The first is based 
directly on an SU{2) Heisenberg magnet, thejSecond on 
"large-N" generalisations thereof. The formeia uses the 
abovementioned small overlap to generate a perturba- 
tion expansion. The valence bond states can be labeled 
by orthogonalised dimer configurations, the quantum dy- 
namics of which is captured by the leading order dimer 
plaquette resonance move, T. The leadpg order diagonal 
term is given by V. The second pathcl generalises the 
SU{2) - Sp{l) to SU{N) on bipartite lattices or Sp{N) 
generally. The latter method essentially generalizes the 
Schwinger boson representation for the spins by intro- 
ducing a large number of additional boson fiavours. In 
the limit N 00, taken at a fixed number of bosons 
per site, the ground states at leading order can be la- 
belled by dimer configurations. It is the next order, 
terms which generate the abovementioned dimer reso- 
nance move. 



III. THE QUANTUM DIMER MODEL AS AN 
ISING GAUGE THEORY 

In the following, we discuss the relationships between 
a number of different models of interest in the context of 
high-temperature superconductivity and quantum mag- 
netism. The naming conventions for them are depicted in 
Fig. ^ Our first mapping - of the QDM to an Ising gauge 
theory (IGT) - proceeds as follows. The naive Hilbert 
space (inclusive of gauge equivalent states) of any IGT is 
defined by an Ising variable tr = ±1 on each link of the 
lattice; each variable will be taken to be the eigenvalue 
of an operator on the corresponding link. We can 
identify the link variable with the presence or absence of 
a dimer on the link, i.e. the number of dimers on each 
link is now given byn = (l + cr)/2 and the dimer number 
operators are fi = ^{1 + a^), where we have suppressed 
the link index. 

Evidently, the naive Hilbert space is too big and we 
must identify the physical subspace that corresponds to 
allowed hardcore dimer coverings, which is done by im- 
posing a constraint at every site that only one link ema- 
nating from it be occupied by a dimer. This is expressed 



as an operation, G, which leaves invariant only the phys- 
ical states, Iphys), those fulfilling the hardcore condition. 
In terms of the operators ct^ , the hardcore constraint be- 
comes: 

^anphys) = (-n, + 2)|phys), (2) 
+ 

where the sum is over all the links emanating from a 
given site i, and ric is the coordination of that site. This 
implies 

Ga Iphys) = Iphys) (3) 

where 




for any a at each site. 

To write the Hamiltonian at half filling, Hj, in a form 
such that [G, Hj] = 0, we define the usual spin-half rais- 
ing and lowering operators ± ia^ , which respec- 
tively add and remove a dimer on a link. We obtain: 

Hi = -tf + vV = -t'^{a+a^a+c7- + h.c.) (5) 

□ 

□ 

where the sites in the last term are labelled sequentially 
around the plaquette, □. The generalization to other 
lattices follows the prescription for writing down QDMs 
outlined earlier. 

Invariance of Hj under local Ising gauge operations is 
easily checked. In fact a larger, U{1), symmetry arises be- 
cause G~^a^Ga = exp(±ia)(T^ so that the phases picked 
up by the products in T cancel. As we will discuss fur- 
ther below, this local U (1) is a consequence of local dimer 
number conservation. So we have the situation that while 
the Hilbert space is that of an Ising gauge theory, the 
physical states and Hamiltonian are invariant under a 
set of continuous local gauge transformations that have 
the form of a U{1) gauge theory. This is a local gauge 
theory version of the more familiar quantum S = 1/2 
XY-model, which also lives in an Ising Hilbert space. 

When holes are doped into the valence bond manifold, 
we need to worry about the comparison of the (poten- 
tially) large hole hopping energy t to the existing mag- 
netic scales that are responsible for the triplet gap. In the 
spirit of RVB theory, we assume that one can approach 
the problem from small values of t and that it is therefore 
sufficient to keep states in which the spins pair up into 
valence bonds except at the sites from which they have 
been removed. The effect of the hole kinetic energy is to 
move a hole and a dimer in ways that keep within this 
manifold. H 

In the Ising gauge language, we have to add a matter 
(Higgs) part, Hm, to the Hamiltonian. Let the Ising 
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G=+l 



GDM 



evenlGT ^ even QDM 

odd IGT — — ^ (odd) QDM 

potential |lenn 

extended QDM 



FIG. 1. Naming conventions used in this paper. Theo- 
ries without matter are referred to as pure, with matter as 
doped. The even IGT is dual to a ferromagnetic transverse 
field Ising model, the odd IGT to a fully frustrated one (see 
Appendix ^) , with the F ^ oo limit corresponding to a pro- 
jection onto the magnetic ground state(s) of the dual Ising 
models. 



Higgs field be denoted by an operator Tx its eigenvalues 
±1 denoting the presence/ absence of a hole. Since each 
site has either a dimer ending on it or is occupied by a 
hole, the constraint equation is modified to 

G Iphys) = exp (^a{f^ + ^(o-"= + 1 - l/nj)^ |phys) 
= Iphys) , (6) 
and the full Hamiltonian is given by iJ = Hi + Hm, with 



H„ 



-m 



(7) 



(ijk) 



and the sum on (ijk) running over triplets of neighbour- 
ing sites. 



IV. ISING GAUGE THEORIES AS 
GENERALIZED DIMER MODELS 

We now review the converse logic of starting with the 
more familiar action of the standard IGT (referred to as 
even IGT from hereon) and ending up with a Hamiltonian 
formulation. The resulting theory is a generalized dimer 
model (GDM) which features dimers on links that obey 
a generalized dimer constraint - in the even case this 
requires an even number of dimers to emanate from each 
site. We also show that the addition of the Polyakov loop 
term Sp to the standard action introduced by SF leads 
to a Hamiltonian formulation that is a GDM with the 
constraint of an odd number of dimers per site and that 
in a special limit, it reduces to precisely a quantum dimer 
model with the hard core dimer constraint. 



A. Hamiltonian vs. Lagrangian formulation 



made anisotropic by choosing a coupling, for plaque- 
ttes containing links in the imaginary time (temporal) di- 
rection, different from that for purely spatial plaquettes, 
£3 This is neccessary to take the time continuum limit 
needed in the derivation of the Hamiltonian. 

We can now choose a gauge wherein all in the time 
direction are -1-1, so that the first term in Eq. || becomes 
a simple bilinear, —K"^ a^a^ , involving only the links 
in the space directions. (Strictly speaking, it is not pos- 
sible to do this as it would have the effect of modifying 
the gauge invariant products of along temporal loops. 
However this obstruction is not important in the time 
continuum limit in the case of the even IGT at T = 0.) 
One then establishes the equivalence to an appropriate 
Hamiltonian and set of constraints in one dimension less 
by comparing the expressions for the partition function 
generated by this action to that arising from a Trotter- 
Suzuki decomposed path integral formulation generated 
by the Hamiltonian 



DM 



(9) 



where the first sum runs over all links and the Hilbert 
space is limited by constraints (below). One finds an 
equivalence of the partition functions in the limit K'^ — > 
oo, with = Kexp(— 2^^^^) fixed. 

The Hamiltonian defined in Eq. ^ retains a gauge in- 
variance under flipping all spins (in the basis) ema- 
nating from one site. This transformation is generated 
by GiGT = ri-i- To reproduce the physics of the even 
IGT we need to impose Gigt = 1 at every site. 

While the Hamiltonian and Hilbert space are naturally 
derived in the basis, the meaning of these constraints 
becomes transparent by considering the system in the cr^ 
basis. For GigtCO |phys) = + |phys), an even number of 
links emanating from site i has = 1. Identifying, as 



we did in Sect. HI, the presence (absence) of a dimer with 
(T^ = ±1, we see that the constraint Gigt — +i implies 
the presence of an even number of dimers emanating from 
each site - whence our label "even" for the IGT under 
consideration. 



B. The odd IGT and the Polyakov loop term Sp 

To obtain a Hamiltonian problem in which the physical 
states have an odd number of dimers at each site (the 
odd IGT), we need to add the Polyakov loop term to the 
action above: 



(10) 



Consider the action of the standard IGT, hereafter re- 
ferred to as the even IGT, 

S = -X^ J2 ^"^"'^'<^' - E f^^f^^^'cr^ • (8) 



where the product runs over all temporal links. This is 
equivalent to assigning a sign to each space-time |-cpn- 
figuration which is the product of Polyakov loopsE3 in 
the temporal direction that wrap around the system for 
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(a) 




H 

















(c) 



FIG. 2. (a) Top: Allowed hardcore dimer configuration. 
Bottom: Allowed configuration in an odd IGT but not QDM. 
(b) Hopping process of hole (denoted by square) generated by 
Hu, which does not conserve dimer number. Two applications 
of Hu yield an allowed final QDM configuration (c); the net 
hop is generated directly by Hm- 



each spatial site. It can be shown (see Appendix ^ that 
this is equivalent, in the time continuum limit, to choos- 
ing GiGT = — 1 in picking physical states for the ac- 
tion of the Hamiltonian Eq. ^. We should note, that 
even with isotropic couplings in the action Eq. ^, Sp 
breaks the symmetry between space and time (lattice 
Euclidean invariance). Consequently one may need to 
be careful about distinguishing the behavior of Wilson 
loops in space and those in time (Polyakov loops) in dis- 
tinguishing confined and deconfined phases - the latter 
are then the correct quantity to calculate. 

For the square lattice, the inclusion of Sp (which arises 
in the work by SF for a Mott insulator with an odd 
number of electrons per site) thus represents a mixtuta 
of dimers (one link occupied) and tetramers (three). EHl 
Whereas the dimers are amenable to an obvious physical 
interpretation as valence bonds, we are not aware of any 
similar interpretation of more complicated polymers. 

We thus see that the somewhat unconventional form 
of the kinetic term in Eq. ||, which consists of raising 
and lowering operators, rather than simply J^ct^, arises 
from the desire to preserve the hardcore dimer constraint 
manifestly. 



to the Lij) appear in the kinetic term T . 

In the presence of doping, there is an additional differ- 
ence between the two theories which lies in the nature of 
the allowed hopping terms. The GDM admits terms of 
the form 



h.c. 



(11) 



corresponding to the processes depicted in Fig. ||b. (The 
sum J2 runs over all sites of the lattice). Note that 
a notion of charge conservation and hence a global U{1) 
invariance, mandates the use of f ^ operators rather than 
f^. In the r — > CX3 limit, this term becomes ineffective 
since one of the two configurations which connects is 
projected out. However, for T large but not infinite we 
generate the term Hm (Eq. ^ with m u'^/T to obtain 
dynamic holes (see Fig. |^). 

Quite generally, extended QDMs thus arise as liimits 
of odd IGTs with additional couplings. /^From the per- 
spective of IGTs, this is a simplification which focuses 
attention on the existence of larger local U{1) invariance, 
but not much more. However from the perspective of the 
physics of antiferromagnets the QDM along with its as- 
sociated conservation of the number of dimers (valence 
bonds) per site appears to be the more natural construc- 
tion. It also has the appeal that it allows two question 
basic to setting up a gauge theory to be answered trans- 
parently: 1) what is the link variable? and 2) what is the 
local constraint? As we have noted, in the QDM limit 
this leads to an Ising link variable which is simply the 
number of valence bonds but a constraint on their num- 
ber more appropriate to aU (1) theory. At the very least, 
short ranged RVB theory is an example of an IGT of a 
strongly correlated system and we will use it to exam- 
ine some of the observations made by SF about IGTs in 
general. More generally, it seems to us that in order to 
make the case for an IGT description in other phases, it 
would be extremely useful to have a comparable identi- 
fication of the link variable and the constraint (even in 
some appropriate coarse grained sense) . 



C. QDM limit of odd IGT 

One can nonetheless retrieve the hardcore constraint 
by explicitly removing the supernumerary dimers 'by 
hand', through a very large coupling constant F. In this 
limit (the F — > oo limit), the term Yia becomes equiv- 
alent to the kinetic term of Eq. ||; it is in this limit that 
the IGT of Senthil and Fisher including the Sp term is 
equivalent to our QDM. In this context, it is interesting 
to note that ite original U(X) gauge theory of Fradkin 
and Kivelsont§ is close in spirit to the above construc- 
tion. There, the presence of a dimer is encoded by an 
angular momentum variable on each link, L^j, which is 
restricted by an analog of the transverse field term to val- 
ues or 1. The raising and lowering operators (conjugate 



D. Duality to Ising models in d = 2 -I- 1 

It is weJl .known that IGTs in rf = 3 are dual to Ising 
models.ll3ll3 In the Hamiltonian formulation in d = 2 + 1, 
there is a simple, geometrical, way of seeing this which 
demonstrates how the standard constraint G = -1-1 trans- 
lates into a dual ferromagnetic transverse field Ising 
model while the alternative constraint G = — 1 translates 
into a fully frustrated Ising model in a transverse field. 
An account of this mapping is given in Appendix ^ 

Evidently, it is most economical to study these dual 
Ising models. An extensive study of frustrated Ising mod- 
els in d = 2 + 1 is reported in Ref. |9^. 
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V. PHASE STRUCTURE AND QUASIPARTICLE 
FRACTIONALIZATION 

For all the pure (undoped) QDMs, the question of pri- 
mary interest is whether they possess a dimer liquid or 
RVB phase. Such a phase automatically leads to free 
spin-1/2 excitations (spinous) and to the decay of an 
ejected electron into a spinon and a spinless charged hole 
(holon) which provides an example of spin-charge sepa- 
ration (in an insulator) in general dimensions. 

The physical arguments leading to the above conclu- 
sions are simple in the valence bond language. A va- 
lence bond can be broken up into two neighboring spins 
1/2. In a valence bond liquid the cost of separating these 
two objects to infinite separation will be finite - hence 
the existence of a spinon spectrum above the triplet gap. 
Further, at large separations one can remove one of the 
spinous to obtain a spinon-holon pair that has the quan- 
tum numbers of a missing electron, or hole. Hence in 
a photoemission experiment one will see a fractionalized 
spectrum above the charge gap. 

Within the framework of the dimer models both 
spinons and holons are represented by monomers and the 
issue is one of computing the free energy of the system 
as a function of monomer seperation. In a liquid phase 
this will be finite. 

It is worth digressing a bit and noting the translation 
between the standard gauge theory lore and the above 
statements. The standard diagnostic of confinement in a 
pure gauge theory is the Wilson loop. In a Lorentz invari- 
ant theory its orientation does not matter and hence we 
may compute the expectation value of a spatial loop as 
well as a temporal loop which is directly related to the en- 
ergy of two separated quarks. In dimer models one does 
not have Lorentz invariance - trivially for we are in the 
time, but not space, continuum limit and less trivially 
because the Polyakov loop term is not Lorentz invari- 
ant, as already noted. Consequently one should really 
compute the temporal Wilson loop (the Polyakov loop). 
Nevertheless, as the following calculation shows, liquidity 
implies a perimeter law for the spatial Wilson loop, and 
as explained above, liquidity in the dimer models signals 
deconfinement . 

A second caveat is that in the dimer model it is the two 
monomer free energy that has a clear meaning. As this 
is a state with physical charge 2 relative to the ground 
state, its free energy cannot be computed as a neutral 
vacuum correlator (unlike a quark/anti-quark potential). 
The strict analog of the quark/anti-quark potential is 
the interaction between a monomer and a site with two 
valence bonds (Fig. ^). Presumably the long distance 
interactions in the two cases will track one another. 

With these comments, consider the spatial Wilson loop 
in the dimer limit. The product 

La 

Trw = Y[<j' (12) 

i=l 




FIG. 3. Top: Neutral 'quark-antiquark' pair. Bottom: 
Charge-2 pair. A square denotes a hole. 
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FIG. 4. Evaluation of a Wilson loop, taken around the cir- 
cumference of the displayed region. The string of operators 
nw annihilates the state unless it encounters an alternating 
sequence of occupied and empty bonds (left panel). The links 
occupied by dimers after the action of nw are denoted by 
empty rectangles. Flipping only one plaquette (marked by a 
plus, right panel) leads to the configuration being annihilated 
by nw- 

reduces to the strings 

■Kw = (CTi^cr^--CrJ,_iO-Z, + ^-C-) (13) 

of dimer creation and destruction operators. In tak- 
ing the expection value of ttw, we select pieces of the 
ground state wavefunction that contain precisely the 
dimer strings in Fig. ^ along the selected loop. To esti- 
mate this fraction we appeal to the extensive entropy of 
dimer configurations and to a healing length ^ in a dimer 
liquid. In the liquid we therefore obtain an estimate 

nw e-^(«)^= , (14) 

that exhibits a perimeter law consistent with the lore 
for a deconfined phase Jl3 Here, c(^) is some numerical 
constant depending on the correlation length. Strictly 
speaking wc have carried out an estimate for a ground 
state wavefunction spread equally over all dimer config- 
urations - such states arise at the so-called Rokhsar- 
Kivelson points (see below) - but this will be qualitia- 
tively correct everywhere in a liquid phase. 

We should note that the above considerations are for 
static matter interacting via fluctuating gauge fields or 
dimers. While it is entirely reasonable that a deconfined 
phase at zero doping will continue into a deconfined phase 
at finite doping, the finite doping problem is logically 
distinct and needs to be treated carefully on its own. We 
note, for instance, the well-known result that dynamic 
matter does not permit an area law phase for the Wilson 
loop, even at arbitrarily weak^oupling on account of the 
screening of the gauge force.LlI 

In the remainder of this section, we discuss the prop- 
erties of dimer models on a number of lattices and for 
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varying dimensionality. Besides presenting a number of 
new results, especially in d = 1, we collate several results 
from the literature, several of which are in the guise of 
stacked magnets or transverse field Ising magnets, which 
therefore need translating. We will also refer to results 
on the GDMs or IGTs outside the QDM limit in places. 

A. d = 1 

While d = 1 is special, it is instructive in that it does 
provide an example of fractionalization that though dis- 
tinct from the higher dimensional versions, fits nicely into 
the QDM description. This point was overlooked by SF 
in their analysis of IGTs. 

Consider first the pure even IGT. In this case there are 
only two states in the Hilbert space, those with — 1 
for = — 1 on all links. The Hamiltonian can only count 
the number of dimers, as there is no local resonance move: 

i? = rX:a% (15) 

whence the cr^ = — 1 state is always the ground state. 
Consider introducing two holes at a separation R. The 
constraint now requires that the links between the holes 
carry = 1 which leads to an energy cost linear in R 
and hence confinement. This is the well-known result on 
the purely confining character of the d = 1 (even) IGT. 

Interestingly, the odd IGT behaves very differently! 
There are still two states in its Hilbert space, but they 
consist of states with alternating values of ct^ = — 1 and 
— 1, i.e. dimers and no dimers - evidently there are 
two such states related by a translation. These states are 
degenerate in energy. Consequently the introduction of a 
hole still produces a domain wall between the two phases, 
but two domain walls do not attract - the charge carriers 
are deconfined solitons. This is, of course, the familiar 
mechanism of solitonie, fractionalization from studies of 
conducting polymers.123 What is interesting is that the 
odd IGT captures this mechanism in d = 1 automati- 
cally. 

This is a good place to give a trivial example of the 
difference between microscopic rewritings and effective 
gauge theories. Consider the one-dimensional Heisenberg 
chain with first Ji and second neighbor J2 antiferromag- 
netic interactions. This can be rewritten as an IGT with 
two Ising gauge fields - one for each bond. As is well 
known, for sufficiently large J2 the chain is in the dimer- 
ized (Majumdar-Ghosh) phase where the effective theory 
clearly involves just one Ising gauge field. 

B. d = 2 

Two dimensions is the case of maximum interest in the 
context of theories of cuprate superconductivity. 



Square lattice: A general analysis of the even IGT cou- 
pled to Ising matter has been given a while back by Frad- 
kin and ShenkerBI They showed that two phases exist: a 
deconfined phase with free charges in the spectrum and 
a confincd/Higgs phase. The phase diagram of the odd 
theory is not known in as much detail. What is known 
is that the undoped odd IGT has a confinement transi- 
tion accompanied by translational symmetry breaking as 
the QDM limit is approached. This result follpsvs from 
analyses of the dual transverse field Ising modellli-as well 
as from a map from the QDM to a height modelS Con- 
sequently, the purely kinetic QDM on the square lattice 
gives rise to a valence bond crystal with confined spinous. 

The extended QDM (Eq. |l]) has been studied in detail 
in and found to be ordered for all values of v, except for 
a transition at the Rokhsar-Kivelson (RK) point v = t 
between a staggered {v > t) and a four-fold degenerate 
state {v < t) which is likely a plaquette state close to the 
RK point and then gives way by a firs±[er|dec-tj:ansition to 
a columnar state at large negative ?;.l3'Ea'E3LZl At the RK 
point, the ground state is an equal amplitude superposi- 
tion of all dimer configurations. Spinons are deconfined, 
precisely at the transition only, and confined elsewhere.c2l 
The unusual feature that a critical point intervenes be- 
tween two crystals, finds an elegant explanation in terms 
of height representations: the effective action in the prox- 
imity of the RK point, has the form conjectured most 
completely by HenleyO 

5 - j d^x dr [(drh)^ + p2(V/i)^ + P4{'^^hf] (16) 

with ~p2 on {v/t) + 1 changing sign precisely at the RK 
point. This action accounts for the crystal for v < t which 
is a flat state of the height variable, the critical correla- 
tions and resonon spectrum cj ~ fc^ at the RK point, and 
the staggered state ioi v > t which corresponds to the 
maximum tilt of the height variable. It also accounts for 
two other non-trivial features of the RK point, namely 
that it has degenerate ground states in all winding num- 
ber sectors (the RK point action is insensitive to gradi- 
ents of heights) and that its equal time height correlations 
are logarithmic and precisely those of the classical dimer 
problem - which follows from the observation, 



J Ll!^ + p^k^ \fPik'^ 

Less is known about the doped QDM, or the odd IGT 
coupled to charged matter. A plausible scenario, based 
on a different large N limit than the one that gives ris€| 
to the QDM, has been advanced by Sachdev and Vojtanil 
who find a set of striped states followed by a d-wave 
superconductor. A direct treatment of the doped QDM 
has not been carried out exceptJ^or an early mean-field 
theory by Fradkin and Kivelson,t3 which led to an s-wave 
superconducting saddle point, but which treats electron 
hopping terms which lack crucial phase factors stemming 
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from the microscopies. Further work on the doped QDM 
would certainly be desirable. 

Triangular lattice: The standard IGT on the trian- 
gular lattice (defined by taking products of around 
elementary triangles) is dual to the ferromagnetijC-trans- 
verse field Ising model on the honeycomb latticed Con- 
sequently, it has a phase transition identical to that of 
the square lattice IGT, between a confined and decon- 
fined phase. While we do not know of a detailed analysis 
of the problem with Ising matter coupled to the Ising 
gauge field, we expect that the Fradkin-Shenker analysis 
applies. 

The odd IGT is dual to the fully frustrated Ising model 
on the honeycomb lattice (see Appendix which has 
been studied by Chandra and ouselveslij with some evi- 
dence for weak ordering involving the breaking of trans- 
lational symmetry in a confining phase. It is also possi- 
ble that the confining phase is absent altogther. Indeed, 
in the extended QDM we have shown that there is def- 
initely a liquid phase for a finite range of parameters, 
2/3 ^ vjt < 1, in which spinons are deconfinedo This 
is the only known example of a deconfined phase in an 
IGT in the QDM limit. We note that a recent neutron 
scattering experiment on a triangular magnet, albeit a 
spatially anisotromc one, appears to have detected de- 
confined spinons H 

Other lattices: The results for the QDM can be gener- 
alised to other lattices. The behaviour of quantum dimer 
models on bipartite lattices follows that of the square lat- 
tice. This is a consequence of the equivalence of the clas- 
sical dimer models to a height model,L3H which in c? = 2 
implies critical correlations which result in an ordering 
transition when quantum fluctuations are switched on.t2l 
In particular, this class, includes the hexagonal lattice 
quantum dimer modelMj 

For QDMs on non-bipartite models, no such general 
result is known. Typically, one expects QDMs to be dis- 
ordered and gapped at the RK point v = t, and we sus- 
pect that this will result in an extended disordered phase 
for V t with spin charge separation in analogy to the 
triangular lattice case. For some interesting results on 
depleted lattices see Ref. 

Charge-2 Higgs scalars: The difference between the 
results for bipartite lattice and the triangular lattice can 
be rationalized by a mechanism that has been invoked 
several times in previous work on the subject. In the 
QDM limit, we obtain a theory with a [/(I) gauge sym- 
metry as we noted earlier, albeit one that lives in an IGT 
Hilbert space. |— . 

A celebrated result of PolyakoveS showed that the 
pure compact U{1) gauge theory confines at all values 
of its coupling in 2-f 1 dimensions. This theory, defined 
by the Maxwell Lagrangian in the continuum, is natu- 
rally formulated on the bipartite square lattice - the pla- 
quette product naturally translates into the forraer La- 
grangian. It was argued by Fradkin and Kivelsor£3 that 
the Polyakov argument goes through for the IGT with 
the Sp term. It follows then that we should expect the 



QDM hmit of the odd IGT to be confining. Evidently, 
such an argument does not by itself rule out a deconfined 
phase in the extended QDM and one needs the more de- 
tailed height representation analysis to show that. 

A second important result on U{1) gauge theories is 
due to Fradkin and Shcnker whoj-atudied their phase 
structure when coupled to matter.til They showed that 
a coupling to charge- 1 matter fields does not allow a de- 
confined phase to exist but that a coupling to charge-2 
scalars did allow a deconfined phase to exist. This re- 
sult provides a rationalization of the triangular lattice 
results: to make contact with the standard results on the 
square lattice, one must treat the triangular lattice as 
a square lattice with additional diagonal bonds. These 
additional bonds imply the presence of a charge-2 scalar 
field coupled to the gauge field, which then opens up the 
possibility of a deconfined phase. 

The general attractiveness of the Fradkin-Shenker re- 
sult to workers on two dimensional quantum magnetism 
is evident. Heisenberg models are easily reformulated as 
U{1) gauge theories but in the search for spin liquids 
with deconfined spinons this is an embarassment. Con- 
sequently it has been suggested by many workers; Read 
and Sachdev,E3 as well as Mudry and FradkinE^) that 
the condensation of an appropriate charge-2 scalar field 
would allow a spin liquid to exist. In Jus work on the 
Sp{N) analysis of the triangular lattice,E2l which exhibits 
a disordered phase, Sachdev again argued that fiuctua- 
tions about the saddle point solution had the structure 
of a U{1) gauge field coupled to a charge-2 scalar, for 
essentially the same reason we invoked above. 

C. d > 2 

In d > 2, little is known about the properties of dimer 
models, even in the classical case. Formally, the Pfaffian 
methods used in d = 2 to gain information about the 
classical models break down due to the overwhelming in- 
crease in thj£.number of terms to be evaluated in higher 
dimensions .Lj 

Within the framework of dimer models, it is likely that 
spin-charge separation becomes prevalent because dimer 
models should become more disordered in high dimen- 
sions. This can be rationalised as the hardcore constraint 
cannot be spread out but number of possible orientations 
increases. 

However, the usefulness of quantum dimer models 
for describing the physics of Mott insulators/Heisenberg 
models decreases in higher dimensions as it becomes in- 
creasingly hard energetically to stabilise valence-bond 
dominated configurations against the Neel state. 

Senthil and Fisher have proposed an experiment to test 
their ideas on a fractionalized phase in the cuprate phase 
diagram and have explicitly linked this experiment to the 
notion of topological order invoked previously in studies 
of the fractional quantum Hall effect by Wen and co- 
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workers.ll3 We examine these assertions in the context of 
dimer models, in reverse order. In particular we will be 
interested in the connection between their idjeas-|aH.d the 
topological analysis of valence bond states .CjIj'Oo 

D. Topological order in the quantum Hall effect 

Quantum Hall states do not break any obvious sym- 
metry captured by a local operator; this excludes cases, 
such as quantum Hall ferromagnets, where a symmetry 
may be broken in addition. There are two alternative 
approaches to characterizing quantum Hall states which 
can both be derived from the rewriting of the electron 
dynamics in terms of bosons coupleAtp-pne or more fluc- 
tuating Chern-Simons gauge fieldsE3il3 In the first, one 
focuses on the bosons and characterizes their condensa- 
tion via an infinite particle electron operator which works 
everywhere in the quantum Hall phase that grows out p£ 
the ideal quantum Hall state in a disordered system.EHI 
In the second approach, one integrates out the bosons 
to obtain a purely gauge action, which then contains the 
Chern-Simons term as its leading piece. 

The Chern-Simons term is topological, i.e. it is insen- 
sitive to the metric of the manifold it is defined on. The 
pure Chern-Simons theory, which describes the strict in- 
frared behavior of the quantum Hall system, has a finite 
dimensional Hilbert space with a set of degenerate states 
whose number depends on the topology of the manifold.Ell 

This leads to the notion of topological order - the idea 
that a quantum Hall state can be characterized by its 
"response" to the topology of the underlying manifold. 
Operationally, one imagines computing the exact spec- 
trum in finite volumes and looking for a low-lying cluster 
of states clearly (i.e. parametrically in system size) sep- 
arated from all other states. This works perfectly for 
clean quantum Hall systems on the torus - e.g. there 
are q exactly degenerate states at filling factor 1/q. In 
this case it is also possible to contruct operators, corre- 
sponding to the adiabatic insertion of one quantum of 
flux through the holes that have the effect to transform- 
ing one ground state into another. As these operators 
commute with the Hamiltonian, their failure to leave the 
ground state invariant was interpreted by Wen and Niu 
as the breaking of a topological symmetry (the symme- 
try algebra itself being dependent on the topology of the 
manifold) . 

To summarize: Topological order in clean quantum 
Hall systems at the ideal fllling factors involves, 

a) a ground state multiplet, separated from other states 
by an amount parametrically larger than the splitting 
between them, and with a degeneracy that increases with 
the genus, g, of the manifold as 

b) a topological symmetry algebra containing operators 
that move the system between different members of the 
ground state multiplet 

c) a long wavelength action (the Chern-Simons action) 



that deflnes a theory with a flnite dimensional Hilbert 
space with the same degeneracies. 

In the clean system, quantum Hall states compete with 
the Wigner crystal, or with various charge density wave 
states. The latter pair of states will lead to a higher 
ground state degeneracy, indeed an inflnite degeneracy 
in the infinite volume limit corresponding to the various 
translations of the crystals as a whole. The same is true 
for the quasiparticle Wigner crystals that will form in the 
close proximity of quantum Hall fillings. In such cases 
of broken symmetry Wen and Niu have argued that the 
splitting between different states will be exponentially 
small in the area of the system (the number of moves it 
takes to convert one ground state into another) instead 
of in the linear dimension, as would be expected from 
tunneling processes involving quasiparticles that would 
move the system between different quantum Hall ground 
states. Hence it would take either a direct examination of 
the states or a study of the magnitude of the splitting to 
decide whether the ground state cluster is due to topolog- 
ical ordering or merely a broken translational symmetry 
(on a manifold of fixed genus, such as the torus which is 
what one is likely to study in practice). An alternative 
approach would be to explicitly lift any degeneracies due 
to broken symmetries by the action of small fields. Any 
residual degeneracy would then be topological in origin. 
For instance the application of a commensurate periodic 
potential would lead to the selection of a unique state 
in the Wigner crystal phase while it would reveal the 
underlying degeneracy of the quantum Hall state in the 
case of a quasiparticle Wigner crystal. (We note such a 
procedure would seem logically necessary for the topo- 
logical degeneracy to track the off-diagonal long range 
order which survives in the quasiparticle Wigner crys- 
tal if its phonons are stiff enough.) Similarly in SU{2) 
quantum Hall ferromagnets the introduction of a Zeeman 
term would be required. Our discussion here illustrates 
two more aspects of topological ordering in quantum Hall 
systems 

d) an exponentially small splitting with linear dimension 
can be attributed to the presence of fractionalized quasi- 
particles that can tunnel across a loop and recombine to 
move the system to a new ground state 

e) it is necessary to break all standard, additional, broken 
symmetries explicitly to reveal the underlying topological 
degeneracy. 

Perhaps the "cleanest" as well as the most realistic way 
to single out the topological degeneracy is to include the 
effect of disorder and thereby exaraiiie a quantum Hall 
phase of finite extent. Wen and NiuEZi have offered argu- 
ments that the inclusion of disorder splits the degeneracy 
by an amount that is 0{e~^^^) in the linear dimension L 
of the system, ^ being a disorder correlation length. Their 
analysis, which holds exactly at z/ = \/q neglects any cre- 
ation of quasiparticles in the ground state itself, i.e. the 
quasiparticle spectrum is assumed to remain gapped in 
the presence of disorder. In general this will not be true, 
and certainly away from u = 1/q there will be localized 
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quasiparticles in the ground state that will give rise to a 
gapless spectrum. In a finite volume, this spectrum will 
acquire a gap that is at worst polynomially small in L 
and so if the exponentially small splitting of the ground 
states remains, they should not prevent an identification 
of the ground state cluster. We do not know of any de- 
tailed examination of whether the ground state splitting 
continues to be exponential in this limit - it would ap- 
pear that one cannot merely argue by continuity from 
the gapped case due to the singular closing of the gap en 
route. Neither is it clear that the operators that move 
us between states in the clean case will continue to work 
with randomly localized (or even crystallized) quasipar- 
ticles - here again we do not know if a generalization is 
possible. Finally, we note that insulating states in the 
disordered system are expected to exhibit unique ground 
states that are separated from excited states by, at best, 
polynomially small gaps coming from localized electrons. 

E. Topological Order in IGTs? 

The general idea of SF is as follows (Ref . Pd,c) . The de- 
confined phase posseses Ising vortex excitations (visons) 
that cost a finite amount of energy. As in any gauge 
theory where such excitations are possible, in a multi- 
ply connected geometry these can be placed so that their 
cores inhabit the holes and we can expect these configu- 
rations to be long lived, and in an appropriate order of 
limits they should be truly meta-stable (that is to say, 
infinitely long lived local minima). For the purposes of 
the experiment proposed by SF (see below) this is suf- 
ficient. To make contact with the notion of topological 
order, SF wish to relate the presence of visons threading 
holes to an infinite volume limit ground state degeneracy 
of 2*^ on a manifold with h holes, that can be interpreted 
as the breaking of a topological symmetry. In the follow- 
ing we explore this idea in some detail with cylinders and 
tori as the manifolds of interest - going beyond those in 
genus while retaining a lattice is tricky, especially when 
the gauge theory arises as an effective theory and so we 
will not venture that far afield. We begin with pure gauge 
theories. 

Gauge fields alone: As noted by SF, the even IGT on 
the cylinder at the point F = exemplifies their ideas. 
There are two exactly degenerate states, which can be 
written in the cr^ representation if one does not worry 
about the constraint. These states have two features of 
note: a) that they exhibit a well defined topological flux 
Fz = Ho ^ where the product Ho taken around the 
circumference of the cylinder. takes the values ±1 in 
the no-vison/vison states; and b) that there exists an op- 
erator Fx = n= where the product n= is taken along 
a seam of links, with the seam running along the axis of 
the cylinder (see SF for details). F^ commutes with the 
Hamiltonian and converts one of the states to the other. 
These two operators capture the two ways of looking at 



the degeneracy, either as a consequence of Ising flux or 
that of breaking a "topological symmetry" in which a 
global operator ceases to annihilate the vacuum. At is- 
sue is whether these generalize beyond this special point 
and to Ising gauge fields coupled to matter, especially in 
the QDM hmit. 

Sticking with the even IGT for the moment, we note 
that the degeneracy is exact for T K in perturbation 
theory, for a cylinder of finite width. In contrast, it is 
clear in the opposite limit F ^ i^T that there is a unique 
ground state. This implies that even for a finite width 
cylinder there is a true phase transition en route. We 
note that F^ commutes with H at all values oiV / K and 
that Fz is a natural order parameter for this transition, 
being odd under the action of Fx. Hence the distinction 
between the two phases is indeed captured by the action 
of Fx and by the development of an Ising flux. We should 
note though that Fz is measurable only for finite cylin- 
ders; being a Wilson loop, it goes to zero exponentially 
in (at least) the width of the cylinder at any F 7^ 0. 

These observations should not really surprise for they 
involve a system that is infinite in two space-time direc- 
tions and finite in one and hence are equivalent to those 
concerning the two dimensional IGT at finite tempera- 
ture. Such a theory indeed posesses a phase transition in 
which the Polyakov loop (a Wilson loop taken in the time 
direction) develops an expection value. In the dual rep- 
resentation this is simply the d=2 Ising phase iiansition 
in a d=3 system that is finite in one direction. o 

We return now to the question of working explicitly 
with gauge invariant states, i.e. those that satisfy the 
local constraint exactly. Given a state | ^E*) in the cr^ rep- 
resentation, we can construct a state V \^) that is gauge 
invariant by the action of the projector 

V^\{{l/2)[GiGT{i) + l] 

i 

which commutes with the Hamiltonian. Evidently, all 
gauge invariant observables have the same value before 
and after the projection. While this indicates that our 
earlier description is correct, it hides a subtlety of some 
interest in making contact with earlier work on the topol- 
ogy of RVB states. To uncover this, note that a state 
written explicitly in the basis is automatically gauge 
invariant if it involves only even numbers of dimers at 
each site. All such configurations can be classified by 
winding numbers - one simply asks how many loops of 
dimers cross a fixed line bisecting a set of horizontal 
bonds. For a finite height cylinder, this number is either 
odd and even and the action of the Hamiltonian pre- 
serves this number. Hence the true ground states must 
be purely even or odd. Now the vison and no-vison states, 
when projected, contain both sectors - they correspond 
to taking the linear combinations |even) ± |odd). Hence, 
although they correspond to a different choice of basis in 
the space of the two degenerate states, it is clear that the 
physical choice for the standard Hamiltonian is that of 
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purely even and odd states which were what were invoked 
in earher studies of RVB states. On the other hand, if 
we were to allow Wilson loops of arbitrary length in the 
Hamiltonian (but with exponentially suppressed coeffi- 
cients to preserve effective locality) we would mix these 
states and obtain the vison/no-vison linear combinations 
split by an amount exponentially small in the cylinder 
circumference. In this case what description one would 
take to be the correct topological decomposition in the 
infinite volume limit would appear to be a matter of taste. 

On the torus, there is no true phase transition even 
for the standard IGT. Instead we find an exponentially 
small splitting between four states when the linear di- 
mension L is increased at a fixed coupling corresponding 
to the deconfined phase and a splitting of 0(i°) between 
a unique ground state and the first excited state in the 
confining phase. In terms of the winding number anal- 
ysis, this corresponds to the four different combinations 
of even and odd in either direction. 

We turn now to the case of the pure odd IGT. Here 
it is instructive to work in the tr^ representation. By 
means of the standard device of using a transition graph 
between a given state and a reference statejs one can 
again assign a conserved even/odd winding number to 
each configuration. For odd height cylinders, a horizontal 
translation by one lattice constant interchanges the two 
sectors. Assuming that odd and even height cylinders 
converge to the same infinite height limit, it follows then 
that the ground state must be at least twofold degenerate 
at all T/K for infinite height cylinders. As the V ^ K 
analysis in the representation is identical to that of 
the even IGT except for a different choice of projector, 

V^\{{l/2)[GIGT{^)-l] , (17) 

i 

there is a twofold degeneracy in that region. Unlike in the 
case of the even IGT, there is a large degeneracy in the 
extreme opposite limit, K = Q, where any dimer covering 
of the cylinder is a ground state. For infinite width cylin- 
ders, i.e. in the two-dimensional limit, there is a four-fold 
crystalline degeneracy as noted earlier. How this degen- 
eracy is modified by finite cylinder widths is not clear to 
us at this point. A preliminary analysis of the QDM on 
cylinders indicates that it will exhibit a two-fold degen- 
erate liquid phase that does not break any symmetries 
as well as a two-fold degenerate columnar phase in which 
the columns run along the cylinder axis. Consequently 
at different cylinder widths the K/T — )■ limit may be- 
have differently. We expect that the large circumference 
limit will be characterized by symmetry breaking which 
may either preserve the two-fold degeneracy of the lowest 
lying cluster (the case if the ground state remains liquid 
for all finite widths) or increase it by a further factor 
of two (the case if the ground state becomes columnar 
already or finite widths). In the former case one would 
have to examine the nature of the degenerate states to 
decide what phase they correspond to. 



On the torus the deconfined phase has again a four-fold 
low lying cluster with a splitting of 0(e^^) while the con- 
fining phase will exhibit a cluster of four low lying states 
with a splitting of 0(e~^ ), corresponding to the necces- 
sity of altering the state over its entire volume instead of 
just along a line in the liquid case. (It is worth noting that 
our previous argument about translations implementing 
winding number sector changes implies that there is an 
exact two fold degeneracy due to translational symmetry 
breaking on odd by even tori.) So on the torus one would 
need to examine the size dependence of the splitting or 
the correlations in the ground states to distinguish the 
two four-fold degeneracies from each other. Alternately, 
as in the quantum Hall case one could turn on symme- 
try breaking fields that would lift the degeneracy in the 
crystalline phase but not in the liquid, deconfined, phase. 
We note that in the context of the cuprates, this is the 
case of maximum interest. 

To summarize: The behavior of Ising gauge fields alone 
does display a "family resemblance" to the quantum Hall 
case with regards to points a,b,d and e made earlier. 
With respect to c the fundamentally discrete character 
of this problem makes it unlikely that there is an ana- 
log. That being said we should note that in the QDM 
limit it does not really go beyond the previous analysis of 
RVB wavefunctions in terms of winding number sectors 
- the latter is an analysis in terms of electric fluxes (the 
momenta conjugate to the gauge fields). 

In this regard the really interesting claim of SF is that 
the phase obtained at finite doping is also character- 
ized by topological order. As the even/odd classification 
breaks down upon doping, this would be a feature not 
obtained by the previous analysis. In the language of the 
IGT we must ask what happens when we add matter to 
the problem. 

Gauge fields with matter: We note at the outset that 
this might be expected to differ from the quantum Hall 
case. In the latter the states differ, in a sense, by the 
insertion of integer numbers of flux quanta through the 
holes. By contrast in the IGT problem, the vison will be 
seen by matter fields as half a flux quantum. 

Nevertheless, the effect of the additional flux can be 
exponentially attenuated if the matter fields are gapped 
on their own. The simplest such case is that of the even 
IGT with Ising matter. While no longer commutes 
with _ff , perturbative considerations indicate that in the 
deconfined phase there are two low lying states with a 
splitting that is 0{e~^) at large L, which goes away on 
leaving this phase. So in this case it is indeed possible 
to relate the deconfined phase to a two-fold degeneracy. 
Having identified the two "ground" states, one can test 
them for the presence of flux. With matter present, the 
even and odd sectors are now connected and the states 
will exhibit (small) expectation values of the Wilson loop 
consistent with the presence and absence of a vison. 

One might wonder if it is possible to relate the two 
low lying states by the action of F^- It turns out that 
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the attempt to create one from the other by its action 
will yield a vanishing overlap in the limit of infinitely 
long cylinders. This result can be obtained perturba- 
tively near the trivial point T/K = 0, u/A = 0. At this 
point, the ground state with = 0, |$o) has cr^ = 1 and 
= —1, whereas the state |<i>i) = |$o) differs in that 
the horizontal are flipped along one seam along the 
axis of the cylinder of height H that the lattice resides 
on. Carrying out perturbation theory to second order in 
u/X yields the perturbed wavefunctions |<i>§) and |$i), 
respectively. One then finds 

= cM-2H{u/X)^). (18) 

where we have exponentiated the linear answer that per- 
turbation theory actually produces. General random 
walk arguments indicate that the exponential dependence 
on the height is exact though the coefficient will be mod- 
ified at higher orders in perturbation theory. In sum, the 
degeneracy is recovered in the infinite system size limit 
but the topological symmetry operation no longer takes 
us between ground states. 

The case of greatest interest is that of charged mat- 
ter coupled to Ising gauge fields. SF have suggested that 
spinon and holon fields coupled to an Ising gauge field 
are the correct low energy theory of a variety of strongly 
correlated systems and have argued that anomalous non- 
superconducting phases would be charaterized by topo- 
logical degeneracies that could, in principle, be used to 
search for such phases in numerical studies or variational 
studies. 

In the QDM framework, we are concerned with adding 
holons to a dimer liquid. If the dimers remain liquid, then 
we have a doped phase that might be expected to inherit 
topological degeneracies from the parent insulating state. 
It would appear that there are three possibilities: 

a) the holons localize 

b) the holons are bosonic and condense thereby giving 
rise to a superconductor 

c) the hplpns are fermionic and produce a gapless 
spectrum.Eil 

In case a) one has perhaps the strongest argument for a 
surviving topological classification and associated degen- 
eracy. Certainly if the holons are truly immobile, one can 
define even and odd sectors for that given configuration. 
If they are localize on some length scale, the classifica- 
tion is no longer strict but it seems plausible that for 
system sizes much bigger than their localization length, 
the degeneracy is recovered. 

In case b) the system ends up with a superconducting 
vortex threading it and so the question is moot. 

In case c) we would truly have a non-Fermi liquid but 
metallic phase. Unfortunately in such a system it would 
appear that all gaps are polynomially small and so it 
will not be possible to select a ground state multiplet 
in an operational sense. From the point of view of the 
QDM, all states involve holons and dimers in correlated 
motion around the torus and no topological character is 



evident. As we were unable to construct use the topolog- 
ical symmetry operator of the pure gauge theory iit-the 
case of Ising matter, we will not succeed here either .E3 It 
would appear then that in this case the non- Fermi liquid 
character will not give rise to a meaningful topological 
degeneracy. 

F. Flux trapping experiments 

We are however, still left with the possibility that 
the states of the doped QDM are characterized by fi- 
nite (if exponentially small) Ising flux measured by the 
Wilson loop. If such a state has a net vison content in 
a non-superconducting phase, it would seem likely that 
it would nucleate a vortex if the parameters are changed 
to condepae the holons. This would then realize the SF 
scenario .Ej 

From our considerations in the last section, we con- 
clude that a flux trapping experiment that cycles between 
phases with the holons localized and then superconduct- 
ing would be the most robust while that between the 
latter and the strange metal is hard to predict without 
a more detailed theory of the metal. In either case, the 
issue appears quite delicate from the QDM viewpoint, 
in which the system is required to remember rather deli- 
cate phase relationships between different components as 
the parameters change. Of course one of the strengths of 
the vison viewpoint is, that by focussing attention on the 
relevant collective co-ordinate, it suggests that this is an 
artefact of looking too microsocpically. Further studies 
of the doped dimer model could be very instructive in 
this regard. 

VI. DISCUSSION 

In this paper we have established and discussed several 
important connections existing between short-range RVB 
phases, quantum dimer models, and Ising gauge theories, 
which have significant implications for the problem of 
spin-charge separation in strongly correlated systems. 

To begin with, we showed that there exists a natu- 
ral physical interpretation of the Hilbert space of RVB 
phases, and that its Ising character follows directly from 
the nature of the states themselves: short-ranged RVB 
states are naturally described in terms of short range 
spin singlets which are either present or absent. Thus, 
from the point of view of the space of states, a descrip- 
tion of the dimer Hilbert space should have a natural 
description in terms of Ising variables living on the links 
of the lattice. As a naive description of this form is seri- 
ously overcomplete, it is clearly necessary to impose con- 
straints at each site which then generate a family of local 
gauge transformations that leave the Hamiltonian invari- 
ant. An Ising constraint would be sensitive only to the 
number of valence bonds modulo two. However, since the 
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number of valence bonds (dimers) is conserved, the effec- 
tive Hamiltonians associated with these states must have 
a natural local conservation law and consequently a local 
U{1) symmetry, instead of the Z2 "natural" symmetry of 
an Ising Hilbert space. Wc further showed that quantum 
dimer models can indeed be realized as (odd) Ising gauge 
theories with additional couplings which project out for- 
bidden configurations of dimers (valence bonds). Thus, 
while the Ising variables provide a natural and economi- 
cal description of the Hilbert space, the native symmetry 
to the physics of short-ranged RVB states is actually U{1) 
and not Z2. 

However, we also found that the phase structure of 
generalized quantum dimer models depends on how the 
local C/(l) symmetry is realized. Superficially, a U{1) 
gauge symmetry may seem to rule out deconfined phases 
since it is quite well known that the vacuum sector of 
U{1) gauge theories are confining in 2 -I- 1 dimensions. 
It turns out that for the case of the gauge theoretic de- 
scription of quantum dimer models the situation is more 
subtle. For instance, on the square lattice the ground 
state turns out to be confining, and thus it is not a spin 
liquid. In contrast, on non-bipartite lattices the situation 
is quite different. Indeed in such cases dimers connect- 
ing sites on the same sublattice give rise to matter fields 
that carry two units of the U{1) gauge charge. In this 
case the deconfinement mechanism of Ref. pl^ (derived 
for the even IGT) can be expected to apply and both 
a confining and a deconfined phase may exist. In the 
deconfined phase the effective remaining "unbroken" lo- 
cal symmetry is reduced to Z2. Thus, this mechanism 
of spin-charge separation relies entirely on the existence 
of a deconfined phase in the Ising gauge theopy. A local 
Z2 symmetry is also central to the work of SF|3 although 
their point of departure is a superconducting state with 
Cooper pairs. We have noted that their starting Hamil- 
tonian has more degrees of freedom than the single band 
t— J type models that we have in mind so their identifica- 
tion of the Ising variable is not as microscopic. However 
valence bonds are sufficiently akin to Cooper pairaa that 
one is tempted to guess that both approaches describe 
the same physics. 

The considerations presented above assume that the 
confinement-deconfinement structure of the phase dia- 
gram of even Ising gauge theories holds also for the odd 
Ising gauge theories. Although this is not rigorously es- 
tablished, there is substantial evidence, including the re- 
sults reported in this paper, that the main difference bew- 
teen even and odd theories is to associate confinement 
with phases in which translation and/or rotational in- 
variance are spontaneously broken, such as valence bond 
crystals and stripe states. In contrast, deconfined phases 
are always liquids. The exception to this is the case of 
d = 1. Here the even IGT, whose ground state is transla- 
tionally invariant, confines at all couplings while the odd 
IGT whose ground state breaks translational symmetry, 
and hence would be expected to be confining by our pre- 
vious remarks, allows test charges to be separated at a 



finite cost in energy. This peculiar feature is, of course, 
the topological mechanism of spin-charge separation in 
d — 1 wherein the charges are accomodated on a pair of 
solitons interpolating between the two ground states. 

A conclusion that emerges from this line of argument, 
is that there is a fundamental difference behind the 
mechanism of spin-charge separation in one-dimensional 
and two-dimensional systems. Indeed, in one dimension 
holons and spinous are actually topological solitons, and 
spin-charge separation is a topological phenomenon, pe- 
culiar to the kinematics of one-dimensional systems. In 
contrast, in two dimensions (and higher) spin-charge sep- 
aration relies on the existence of deconfineinent in the 
sense of liquidity, which is a property of the spectrum of 
states in a particular phase of matter, and as such it does 
not hold in general; deconfinement takes place in some 
cases, such as the triangular lattice which can have a spin 
liquid ground stateJ23 whereas-mnfinement is naturally 
realized on the square lattice.Ha 

The question of the existence of a deconfinement mech- 
anism of gauge theories with dynamical matter at finite 
density has a long history in high energy physics which 
which is rather similar to the quest for a spin-charge sep- 
arated state in condensed matter physics. The difficulties 
of defining order parameters and other tests of confine- 
ment has been a central theme in that field since the 
late seventies. In fact it has long been recognized in that 
field that no such tests can exist in terms of gauge invari- 
ant local operators (such as order parameters) or Wilson 
loops, if the dynamical matter fields carry the fundamen- 
tal gauge charge. A related and important current ques- 
tion is if hadronic matter at finite density is generally and 
smoothly connected to conventional nuclear matter, or if 
a genuine quark-gluon plasma exists as a state of matter 
with unique measurable signatures. This latter phase is 
indeed precisely the equivalent of the spin-charge sepa- 
rated phase discussed here. 

Finally, we have also discussed the question of topo- 
logical degeneracy of the deconfined spin liquid states, 
and their possible detection which we argue is not con- 
tingent upon the former in any precise sense. We have 
discussed in some detail the set of desiderata associated 
with the notion of a topological degeneracy by reviewing 
the case of clean quantum Hall systems at the ideal filling 
fractions. We have discussed the applicability of these to 
disordered quantum Hall systems and then to the case of 
Ising gauge theories. We find that while there is certainly 
a sense in which IGTs in their deconfined phases exhibit 
a finite ground state degeneracy in the thermodynamic 
limit, in general there is no accessible operational test for 
this degeneracy short of a full solution of the spectrum 
of states. In particular we find that the overlap of a test 
state with one vison wrapped around a non-contractible 
loop is orthogonal to any ground state in the thermody- 
namic limit, and therefore it does not connect distinct 
degenerate states. This behavior stands in contrast with 
what happens in ideal quantum Hall fluids and chiral spin 
states, although it may be generic in more realistic cases. 
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in the transition graph any line connecting that spin to 
the reference spin crosses. If this number is even, the 
spins are aligned, otherwise they are antialigned. 

To construct the equivalence between the Hamiltoni- 
ans, we only need to note two facts. Firstly, the presence 
of a satisfied bond gains an energy fc, whereas a frustrated 
bond costs the same amount of energy. Translating this 
into a statement about the absence, presence of a dimer, 
we obtain the identification T — k for the first pair of 
coupling constants. Secondly, note that flipping a spin 
Si implies exchanging all its satisfied bonds for frustrated 
ones and vice versa. This is equivalent to exchanging oc- 
cupied and empty dimer links of the plaquette i, at the 
center of which Si is located. This immediately yields the 
identification of the spin flip effected by the S'^ operator 
with the plaquette term in Hgjjm, together with j = n. 
This completes the demonstration of duality. 



APPENDIX A: DUALITY OF IGTs WITH ISING 
MODELS IN d = 2 + 1 

We show that the GDM with the Hamiltonian given 
by Hg dm (Eq- ^) in c? = 2 + 1 is dual to an Ising model 
with the Hamiltonian: 

G = +l: H+ = -kJ2S!S^+lJ2^^ ^ (Al) 

G=-l: //_ = -^fc,,5f5|+7^^f , (A2) 

where the sums run over the links of the dual lattice 

and the ^ over the sites. The S are Pauli spin operators 
defined on sites of the dual lattice, fc > and \kij\ = k. 
The case of G = — 1 is known as a fully frustrated Ising 
model (FFIM) since each plaquette □ has to have at least 
one frustrated bond: Yiai^ij /^) ~ whereas the case 
G = -1-1 is a ferromagnetic Ising model (FIM). 

The starting point of the duality is the identification 
of a frustrated bond in the Ising model with a dimer in 
the GDM. One can easily convince oneself that each pla- 
quette in the FIM (FFIM) has to have an even (odd) 
number of frustrated bonds, which takes care of the con- 
straint G = +l (-1). 

Gonversely, each dimer state corresponds to a unique 
spin state (up to a global Ising reversal). This can be 
seen by taking a reference spin configuration, for example 
S^ = 1, which corresponds to a reference configurations 
of dimers, namely one without dimers (G — -1-1), or to a 
columnar dimer state (G — —1). Any other dimer config- 
uration can then be used to generate a 'transition graph' 
(see Ref. obtained by superimposing that dimer con- 
figuration with the reference dimer configuration. The 
resulting transition graph contains only closed loops. To 
fix the overall Ising redundancy, an arbitrary reference 
spin is chosen to point up.cJ The orientation of any other 
spin is then obtained by counting the number of dimers 



APPENDIX B: THE POLYAKOV LOOP TERM IN 
THE ACTION 

In order to see the connection between the constraint in 
the Hamiltonian formalisip-|aad the role of the Polyakov 
loop in the path integral,E3C3 it is useful first to recol- 
lect the appropriate construction for the electromagnetic 
gauge field. The Lagrangian density for the free electro- 
magnetic field is 

where 

Ei = doAi + diAo 
Bi ~ €ijkdjAk 

and Jo is a static charge distribution, say 

jo{z) = ^{z -x)- 5{z - y) 

for two static charges aX z — x^y (with charge ±1 respec- 
tively). The path integral in D space-time dimensions 
is 

f i I d'^xClA.j] 

Z[j] = J VA^ e J 

r i f d''x[C[A,0]~Aojo] 

^ VA^ e J 

Thus, 

Z\j] / dxoAo{x,xn) +i / dxoAo{y,xo) 

namely, the expectation value of the product of two 
Polyakov loops. 
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It is easy to show that in the Hamiltonian picture the 
Polyakov IpQps become static sources in the Gauss' Law 
constraint Let us rewrite the path integral by using the 
coherent state representation, which is a integral over 
both the vector potential Ai, the conjugate momenta, 
the electric field Ei, and the Lagrange multiplier field 
Aq. Glossing over issues related to gauge fixing, gauge 
copies, and Faddeev-Popov determinants, one writes. 



Z[j] = / VE,VA,VAo e 



£[Ai,E„Ao,j] 



where 



+Ao {d,E,~jo) 

Thus, we see that the role of Aq is of a Lagrange multi- 
plier that forces Gauss' Law 







as a constraint on the physical Hilbert space. Thus, the 
Polyakov loops are equiv alent to static sources. Notice 
that this is really the Hamiltonian picture since we get 
that the momentum canonically conjugate to A is — -B, 
as we should. 

We now turn to the case of the Ising gauge theory. 
Consider the Hamiltonian of Eq. ^. For convenience, we 
define /i^ = (1 — (t^)/2, so that the Hamiltonian is writ- 
ten, up to a constant. 



□ 



(Bl) 



It turns out to be convenient to rewrite the constraint 
operator G = YI+ '^^ follows. Let L{i) = v{i)'-^ + 
fv^ , where v^{i) = (1 — 'r^(z))/2. Then the projector 

enforcing G{i) |phys) — T(i) |phys) at site i is given by 



Pii) = (1/2) h+(-i)(i'W+€W) 



(B2) 



Here, £,{i) = (1 -I- T{i))/2 is for an even site and —1 
for an odd site. In the absence of matter, all sites of the 
even (odd) theory are even (odd), but the addition of a 
hole at a site changes it from even to odd and vice versa, 
so that the following treatment is appropriate for static 
matter. 

We now Trotterise the partition function Z{(3) at tem- 
perature 1 //? and obtain 



Z{/3) = Tr (^cxp(-/3i7)p) 



lim n {{al^,}\cxpi-eH)P\{al}) , (B3) 
C=o 



where the Greek letter ( labels the (imaginary) time 
shces, and e = P/N, and the tr^ are eigenstates of . 

Consider a single term in the product, which we eval- 
uate by inserting a complete set of eigenstates of a^: 

({a,Vj|exp(-e7?)P|{a|}) =Tr|,.} ^ (1/2)^'- 



{\{z)=Q.l} 



exp 



n 



^-EK^cr^a^aV^^ exp ^-2er^/x=^^ 
!xp ^7^A^(^) 



(K+i}IK})(K}IK}) 

Here, is the number of sites, and Ni, the number of 
bonds. We have rewritten the projector as an expo- 
nential and turned the operators into numbers by let- 
ting them act on their appropriate eigenstates. Note 

that (^{a^} = 2-^>-/^exp(iTTj2-tJ-c^c)' where 

the sum runs over all links in timeslice C. 

Collecting together the terms involving the cr^ , we ob- 
tain 

J2 n ^e'^P {^c(*^') [-2er + iTT (Xdi) + Ac(^ + V) + 

m!M+M?+iM)]} • (B4) 

Here, zp labels the bond connecting site i with its neigh- 
bour in a spatial direction labelled by V. 

The term in parentheses in the foregoing equation can 
be turned into a plaquette product by defining cr^ — 
(1 + A)/2 on the temporal links, so that this expression 
becomes 

i-D \ D i-D 

In the last step, we have used the fact the product over 
plaquettes containing temporal bonds occurring in this 
expression, Yin '^^ i can only take on values ±1. The 
equality holds for a temporal coupling K'^ = tanh(— 2er). 

Putting this result back into Eq. B3, using 2 cosh K'^ 
exp (if '^) for e — > 0, and substituting for A in terms of 
temporal cr^, we obtain 



Z{f3) = (1/2)^-^ Tr 



exp 



□ 



(B6) 



where = ek, the first sum in the second line runs over 
spatial plaquettes, the second over temporal plaquettes. 
The trace now runs over all the cr, both in spatial and 
temporal directions. 
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Crucially the product JJ^ runs over the temporal links - 
this is the Polyakov loop term. It contributes a nontrivial 
phase for all the odd sites. This is what we set out to 
show. 
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